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Abstract. We performed numerical simulations of dynamical encounters between hard massive 
binaries and a very massive star (VMS; formed through runaway mergers of ordinary stars in 
the dense core of a young massive star cluster), in order to explore the hypothesis that this 
dynamical process could be responsible for the origin of high- velocity O 200 - 400 kms -1 ) 
early or late B-type stars. We estimated the typical velocities produced in encounters between 
very tight massive binaries and VMSs (of mass of > 200 M©) and found that about 3 — 4 per 
cent of all encounters produce velocities of ^ 400 kms , while in about 2 per cent of encounters 
the escapers attain velocities exceeding the Milky Ways's escape velocity. We therefore argue 
that the origin of high- velocity (^ 200 — 400 kms -1 ) runaway stars and at least some so-called 
hypervelocity stars could be associated with dynamical encounters between the tightest massive 
binaries and VMSs formed in the cores of star clusters. We also simulated dynamical encounters 
between tight massive binaries and single ordinary 50 — 100 Mq stars. We found that from 1 to 
~ 4 per cent of these encounters can produce runaway stars with velocities of ^ 300 — 400 kms -1 
(typical of the bound population of high-velocity halo B-type stars) and occasionally (in less 
than 1 per cent of encounters) produce hypervelocity (^ 700 kms ) late B-type escapers. 
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1. Introduction 

The origin of high-velocity runaway stars can be attributed to two basic processes: 
(i) disruption of a tight massive binary following the (asymmetric) supernova explosion 
of one of the binary components (Blaauw 1961) and (ii) dynamical three- or four-body 
encounters in dense stellar systems (Poveda, Ruiz & Allen [1967) . In the first process, 
the maximum velocity attained by runaway stars depends on the magnitude of the kick 
imparted to the stellar supernova remnant [either a neutron star (NS) or a black hole 
(BH)] and for reasonable values of this magnitude, the runaway velocity does not exceed 
~ 200 kms -1 (e.g. Portegies Zwart [2000; cf. Gvaramadze 2009). In the second process, 
the ejection velocity could be higher. For example, the maximum velocity that a runaway 
star can attain in binary-binary encounters is equal to the escape velocity from the surface 
of the most massive star in the binaries and could be as large as ~ 1400 kms" 1 (Leonard 
fTMll . 

The recent discovery of the so-called hypervelocity stars (HVSs; e.g. Brown et al. 
2005J - the ordinary stars moving with velocities exceeding the Milky Way's escape 
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velocity - attracted attention to dynamical processes involving the supermassive BH 
in the Galactic Centre (e.g. Gualandris, Portegies Zwart & Sipior 120051 Baumgardt, 
Gualandris & Portegies Zwart |2"0"0T)|) . These processes [originally proposed by Hills (1988) 
and Yu & Tremainc (2003)] can result in ejection velocities of several 1000 km s _1 . Similar 
processes but acting in the cores of young massive star clusters (YMSCs) and involving 
dynamical encounters with intermediate-mass (~ 100 — 1000 M Q ) BHs (IMBHs) were 
considered by Gvaramadze, Gualandris & Portegies Zwart (2008) to explain the origin 
of extremely high-velocity (~ 1000 kms" 1 ) NSs (e.g. Chatterjee et al. [2005]) and HVSs. 
Gualandris & Portegies Zwart : 2007; proposed that exchange encounters between hard 
binaries and an IMBH formed in the core of a YMSC in the Large Magellanic Cloud 
could be responsible for the origin of the HVS HE 0437— 5439. A strong support to the 
possibility that at least some HVSs originate in star clusters rather than in the Galactic 
Centre comes from the proper motion measurements for the HVS HD 271791, which 
constrain the birth place of this star to the outer parts of the Galactic disc (Heber et al. 
12008)1 . 

2. High-velocity runaway stars from three-body encounters 

In this paper, we explore the hypothesis (Gvaramadze I2007|) that some high- velocity 
runaway stars could attain their peculiar velocities in the course of strong dynamical en- 
counters between hard massive binaries and a very massive (^ 100 — 150 M Q ) star (VMS), 
formed in the core of a YMSC through collisions and mergers of ordinary massive stars 
(e.g. Portegies Zwart & McMillan 12002]) . In this process, one of the binary components is 
replaced by the VMS, while the second one is ejected (the so-called exchange encounter), 
sometimes with a high velocity. Our goal is to check whether or not this process can 
produce early B-type stars (the progenitors of NSs) with velocities of 200 — 400 kms" 1 
(typical of pulsars) and 3 — 4 M Q stars with velocities of 300 — 400 kms^ 1 (typical of 
some late B-type halo stars). 

In our study we proceed from the similarity between encounters involving a VMS 
and those involving an IMBH (the latter process is already known to be able to pro- 
duce high-velocity runaways; Gvaramadze et al. 120081) and the fact that the radii of 
VMSs are smaller than the tidal radii of the intruders (the massive binaries), so that 
the tidal breakup and ejection can occur before the binary components merge with the 
VMS (Gvaramadze 2007; Gvaramadze et al. I2009( ) . Our study is motivated by the recent 
evolutionary models of VMSs (Belkus et al. I2007| Yungelson et al. I2008j) . which sug- 
gest that VMSs can lose most of their mass via copious winds and leave behind IMBHs 
with masses of $J 70 M©, which are not large enough to contribute significantly to the 
production of high- velocity runaway stars (see Gvaramadze et al. 2008). We therefore 
explore the possibility that a VMS could produce high-velocity escapers (either early or 
late B-type stars) before it finished its life in a supernova and formed a BH. To check 
this possibility, we performed numerical simulations of three-body exchange encounters 
using the sigma3 package, which is part of the STARLAB software environment (McMillan 
& Hut 119961 Portegies Zwart et al. I2001|) . For a detailed description of the simulations 
see Gvaramadze et al. (2009]). 

First, we focus on exchange encounters producing high- velocity early B-type stars 
and simulate encounters between a binary consisting of two main-sequence stars with 
masses mi = 40 M Q and m-i — 8 M© and a massive compact body, either a VMS or 
an ordinary star of mass of 50 — 100 M© (the most massive ordinary stars formed in 
clusters with mass of ~ 10 3 - 10 4 M Q ; e.g. Weidner, Kroupa & Bonnell 2009). In order 
to maximize the ejection speed, we assume that the binary system is very tight, e.g. 
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Figure 1. Left: Velocity distributions at infinity for escaping ( 8 Mq) stars in encounters between 
a binary consisting of a primary star with mass mi = 40 Mq and a secondary star with mass 
77i2 = 8 Mg, and a single VMS of mass M V ms = 50 M , 100 M Q , 500 M Q and 1000 M (from 
left to right). The binary semi-major axis is a = 0.15 AU. Right: The probability of exchange 
encounters between a (40,8) Mq binary (with a = 0.15 AU) and a single VMS resulting in 
ejection of the 8 Mq star with a velocity from 200 to 700 kms" 1 (top to bottom). 
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Figure 2. Left: Velocity distributions at infinity for escaping stars in encounters between a 
binary consisting of a primary star with mass mi = 10 Mq and a secondary star with mass 
m 2 = 4 M Q , and a single VMS of mass M V ms = 50 M , M V ms = 100 Mq, 500 M Q and 
1000 Mq (from left to right). The binary semi-major axis is a = 0.06 AU. Right: The probability 
of exchange encounters between a (10,4) Mq binary (with a = 0.06 AU) and a single VMS 
resulting in ejection of the 4 Mq star with a velocity from 200 to 700 kms -1 (top to bottom). 



formed by tidal capture. In this case, a ~ 3ri ~ 0.15 AU, where r\ is the radius of 
the primary star (see Gvaramadze et al. I2009|) . Fig.Q] (left panel) shows the velocity 
distribution for 8 M Q escapers. As expected, the more massive VMSs are more likely to 
eject stars with high velocities. In the right panel of Fig.Q] we show the probability of 
exchange encounters resulting in ejection of the 8 Mq binary component with a velocity 
from 200 to 700 kms -1 (top to bottom). For Mvms ^ 100 Mq about 3 per cent of all 
encounters produce runaways with peculiar velocities ^ 400 kms -1 . It can be seen that 
even an ordinary star of mass of 50 Mq can occasionally (in ~ 1 per cent of encounters) 
produce an escape velocity 400 kms -1 . In order to produce escapers with velocities 
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typical of HVSs (Vac ^ 700 kms -1 ), a VMS of several hundred solar masses is required. 
In the case of a 200 — 300 Mq VMS, ^ 2 per cent of all encounters result in an escape 
velocity of ^ 700 kms . This fraction gradually increases to 3 per cent for the more 
massive VMSs. 

Next, we consider exchange encounters producing high- velocity late B-type stars and 
simulated encounters between a very tight binary (mi — 10 M , m 2 = 4 M Q and 
a ~ 3ri ~ 0.06 AU) and a (very) massive star of mass of 50 — 1000 M Q . The velocity 
distributions for 4 M Q escapers are shown in Fig. [2] (left panel) for four different values 
of the VMS mass: M VM s = 50,100,500 and 1000 M . The right panel of Fig.0 shows 
that for all values of Mvms about 2 per cent of encounters results in peculiar velocities 
^ 300 — 400 kms -1 and that about the same percentage of ejected stars attains velocity 
^ 700 kms" 1 if Mvms ^ 200 M Q . 

We therefore argue that the origin of high- velocity (^ 200 — 400 kms -1 ) early and late 
B-type runaway stars and at least some HVSs could be associated with dynamical en- 
counters between the tightest massive binaries and VMSs formed in the cores of YMSCs. 
The future proper motion measurements for HVSs with GAIA will reveal what fraction 
of these extremely high-velocity stars originated in the Galactic disk. 
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